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We consider the problem of the initial conditions and behaviour of the perturbations in scalar field 
cosmology with general potential. 

We use the general definition of adiabatic and isocurvature conditions to set the appropriate initial 
values for the perturbation in the scalar field and in the ordinary matter and radiation components. 
In both the cases of initial adiabaticity and isocurvature, we solve the Einstein and fluid equation 
on : at early times and on superhorizon scales to find the initial behaviour of the relevant quantities. In 

ON . particular, in the isocurvature case, we consider models in which the initial perturbation arises from 

the matter as well as from the scalar field itself, provided that the initial value of the gauge invariant 
. curvature is zero. 

D ' We extend the standard code to include all these cases, and we show some results concerning 

the power spectrum of the cosmic microwave background temperature and polarization anisotropies. 
[ — ■ In particular, it turns out that the acoustic peaks follow opposite behaviours in the adiabatic and 

isocurvature regimes: in the first case their amplitude is higher than in the corresponding pure cold 
dark matter model, while they make the opposite thing for pure isocurvature initial perturbations. 
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I. INTRODUCTION 



In recent times, the need for a "quintessence" component has come out due to the several difficulties of the standard 
f2 m = 1 Cold Dark Matter (CDM) model which proved unable to explain the observed features of large scale structure. 
' In the context of inflationary cosmologies, we expect that the present spatial curvature of the Universe is negligible and 
the total energy density equals the critical energy density; on the other hand, there is growing observational evidence 
that the matter energy density is remarkably below the critical value, even taking into account exotic and so far 
■ undetectable particles known as cold dark matter. Thus, we are faced with figuring out how to explain missing-energy 
Qh values of as much as 70 or 80 percent of critical density. 

Further, there is need to have the age of the Universe, to, exceeding the age of globular clusters in our galaxy; 
the limits on to are holding at about 13 Gyr or more M, and when combined with current estimates of the Hubble 
expansion parameter, converging from different methods to Ho ~ 60 ± 10km/ sec/ Mpc give rise to an observed 
value of the "expansion-age" parameter ifo^o — -8, sensibly higher than 2/3 as predicted by the standard Einstein-De 
Sitter model. 

Preserving the flatness of the Universe, its age could be enhanced by lowering the matter content in models 
involving a component whose equation of state is different from matter and radiation, for example in models including 
a cosmological constant. One more problem with the CDM model arises from the mismatching of the galaxy clustering 
power spectrum shape, when only COBE normalized spectra are considered. All this motivations have ratified the 
demise of the standard CDM model, leaving cosmologists with the open question of what the missing-energy candidate 
could be. 

Apart from the cosmological constant, introduced in some models in order to maintaining the spatial flatness, but 
still retaining serious unsolved theoretical issues, several "quintessence" models where proposed as candidates for 
the missing energy, often modeled as scalar fields rolling down their potentials 0-0], or mo re g enerally described 
in terms of an unspecified equation of state different from that of matter and radiation p2-n4T. We refer to any 
component whose properties are well described in terms of a scalar field evolving in a potential which couples to 
ordinary matter only through gravity. In some sense, it can behave like a cosmological constant when its kinetical 
energy is negligible with respect to the potential energy, so that the scalar field equation of state approaches -1; due 
to the strongly relativistic nature of such component, the characteristic scale of clustering processes for a scalar field 
is just the horizon, [ |ll| giving a similarity with a cosmological constant in the undetectability of quintessence energy 
concentrations on scales smaller than the horizon. 

The interesting feature of the "quintessence" component is just that, contrary to the cosmological constant, it is 
time- varying and spatially inhomogeneous, so that it can develop fluctuations which can be relevant in the perturbation 
growth and can leave a characteristic signature in the CMB and in the large scale structure. Even though many of 
these imprints have been caught in previous works, the issue of initial conditions and scalar field perturbations has 
often been underestimated; in particular, we found a gap regarding the opportunity to impose isocurvature initial 
conditions in a several-components system including a minimally coupled scalar field. Our aim is to give a complete 
prescription for describing adiabatic and isocurvature initial conditions if an additional component is present in the 
form of such scalar field; this can be acquired by giving the set of equations relating all the fluid components needed 
in the two cases. 

In order to do that, we need the background and perturbation equations which we briefly review on section II and 
section III. The original results of our work are presented in sections IV and V, where we work in the formalism 
of the synchronous gauge, by generalizing the work of Ma and Bertschinger p4| and finding the super-horizon-scale 
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behaviour of perturbations at early times, starting from initial zero entropy perturbations (adiabatic case) or initial 
zero curvature perturbations (isocurvature case). We express the needed gauge- invariant quantities, namely entropy 
and curvature perturbations, in terms of synchronous perturbations of baryons, photons, massless neutrinos, cold dark 
matter and a minimally-coupled scalar field, as well as metric perturbations. In section VI the results are translated 
in conformal Newtonian gauge, and in section VII we numerically investigate the growth of entropy and curvature 
perturbations starting from different initial conditions, and we compare them with the corresponding behaviours in 
the standard CDM model. In that computation, the adopted scalar field is associated with an ultralight pseudo 
Nambu-Goldstone boson pq |, with global spontaneous symmetry breaking scale / ~ 10 18 GeV and explicit breaking 
scale M ~ 10~ 3 eV; such a field should be acting at present like an effective cosmological constant and dominating the 
energy density of the Universe. Also we plot, and discuss, the pure adiabatic and pure isocurvature CMB anisotropies 
spectra, again making a comparison with the standard CDM. 



II. EINSTEIN AND CONSERVATION EQUATIONS 



We begin by a brief review of an homogeneous Friedmann Robertson Walker (FRW) cosmology in which there is 
the additional contribution coming from a minimally coupled real scalar field <fi evolving in a potential V((f>). 

We consider only models with total Q = 1 in this paper , and we work in conformal coordinates, so that the line 
element is ds 2 = a 2 [— dr 2 + Sijdx l dx : >} where a is the cosmic scale factor and r is the conformal time. 
The scalar field energy density and pressure, associated to the Lagrangian describing the classical behavior of <fi 

C = -\^—g[g^d^d^ + 2V{4>)], (1) 
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follow from expression of the stress-energy tensor 

^ = ^^4^» + 2 W ' ( 2 ) 



and are given by 



2^ + ^)5 P ^Ya 2 



4? + V(<P) ; P<t> = — ^ 2 - V(<P) , (3) 



where the overdot denotes the derivative with respect to conformal time r. The above quantities evolve according 
to the Friedmann equation , in which we separate the contributes of matter, radiation and scalar field to the total 
energy density, 
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a ' 



n 2 =[-\ = ^-a 2 [p m + Pr + P<p } , (4) 



together with the conservation equations 



2H<j> + a 2 ^- = \^-(a 2 <j>) + a 2 V'(cp)=0 , (5) 
or dr 



p n + nHpn = , (6) 

where Ti. is the conformal expansion rate of the Universe, p n is the energy density contributed by radiation (n = 4) 
or nonrelativistic matter (n = 3) and V — dV/dip. Note that, from (|^), the second order Klein-Gordon equation (|5|) 
is equivalent to the conservation law 

P4, = ~3H(p4, +P4,) . (7) 

Including all the modifications due to the additional scalar field component , we shall carry out a fully relativistic 
treatment of the perturbations of this background , based on the notation of Ma and Bertschinger [|24| . We work in 
the Fourier space and we perform the parametrization of the perturbed quantities in the formalism of the synchronous 
gauge , in which the perturbed line element is ds 2 = a 2 [—dr 2 + (Sij + hij)dx l dx^]. Since we are interested here only 
on scalar-type perturbations, the metric perturbations can be parametrized as 
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hij(x., t) — / d 3 ke 



kikj/i(k, r) + (kikj - -5ij)6T](k, r) 



(8) 



with k = kk and /i denoting the trace of hij . 

Note that the synchronous potentials h and r] in fc-space are related to the gauge-invariant variables of Bardeen 
(1980) and * of Kodama and Sasaky (1984) by the relation 



<&A = * = 



1 



h + 6?y + - ( h + 6?) 



(9) 



which allows us to relate h and ry to the gauge-invariant curvature perturbation £ (see 21 1) 



C = -(ft" 1 * + *)/(l + w) + * 



(10) 



being w = p/p; the above expression will be useful in the following. Now, focus on the equations describing the 
evolution of perturbations involving the various components. 

As it is well known, the scalar field can mimic a cosmological constant if its kinetic energy is negligible with respect 
to the potential one. However, a substantial difference is that it admits perturbations around the homogeneous 
solution of (||); in linear theory, they are described by small fluctuations 6<j) and 8<j> around the background values, 
driven by the equation of motion 



\7 2 6cf> + a 2 



,d 2 V 



> + ±4>h = o 



(ii) 



The density, pressure and velocity perturbations for the scalar field are described as usual by the following quantities 

<pS<p 



- -fiT v - 
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& P4> = ^ 5T i 

(P4> +P<I>)6<I> = = ka~ 2 
= , 



V'5(j> 
V'86 



(12) 

(13) 

(14) 
(15) 



<$0 = 

therefore, we define the differential ratio 



= k—r- 

a 2 P<t>+P<i> cj) 



a 2 V 



+ a 2 V't 



— c <b + T 1 4> J 



that differs from the scalar field sound velocity 



by the term 1^ describing the entropy contribution 
It is useful to describe radiation in terms of the coel 



and polarization brightness functions, Ay(k, n, r) and Ap(k, n, r) 

oo 

A T (k, n, r) = 2(-i , )(2I + l)A T/ (fc, r)F ; (k ■ n) 



(16) 



(17) 



(18) 



19|. 

cients characterizing the Legendre expansion of the temperature 



(19) 



1=0 



A P (k, n, t) = J2(- il )( 21 + 1)Api(*, r)P z (k ■ n) 



(20) 
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Their evolution is completely determined by the Boltzmann equations; denoting by ot the Thomson scattering cross 
section and by n e the electron density, we have for photons: 

(21) 
(22) 
(23) 
(24) 

an, a,. ^(A T2 ( 7 ) + A P0 ( 7 ) + A P2 ( 7 )) I Sio + ^ I - A Pi(7 ) , (25) 



$1 
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2ct 7 




A T ;( 7 ) 


fc 

21 + 1 


Api( 7) 


fc 

2Z 1 


where 





-fcA 



7 3 



4 • 8 
— h + - 
15 5 



9 1 

r) - -an e a T aj + —an e a T (A P0 ^) + A P2 ( 7 )) 



-[/A T(Z _ 1)(7 ) - (I + l)A T(i+1)(7) ] - an e cr T A T , (7 ) (/ > 3) 



-[ZAp (i _x)( 7 ) -(1+ l)Ap( i+ i)( 7 )] + 
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5 1 = A T0 , 8 1 = -kA T i , cr 7 = -A T 2 ■ 
The perturbed stress-energy tensor for radiation contributes with the following non-zero quantities: 



(kjkj 



<5T* 



-p 7 <5 7 + £} 



-p 7 CT 7 



(26) 

(27) 
(28) 

(29) 

(30) 



The expansion (|19|) also applies for massless neutrinos; their evolution equations in the synchronous gauge are given 
by the following system: 



-5 V - o v \ , 



1 



8 n 3 , A 4,8. 
2er„ = —Q v - -rKA T3(i ,) + —h + -r] 
15 5 15 5 

k 



) = ^T[[ lA T(i-i)(») - (i + l)A T(m)(l/) ] (I > 3) . 



(31) 
(32) 

(33) 
(34) 



Pressureless Cold Dark Matter interacts only gravitationally with other particles and in the synchronous gauge its 
peculiar velocity is zero; setting 8 C = 0, the evolution of CDM density perturbations is given by 



$c = — 7:h 1 
2 

and the non-zero component of its perturbed stress-energy tensor is 

8Tq = ~p c 5 c . 

Taking into account the coupling between photons and baryons by Thomson scattering, 

5b = -Ob - -^h , 

9b = ~h + c 2 s k 2 S b + ^ L an e o T (Q 1 - Ob) , 
a 3p b 



(35) 

(36) 

(37) 
(38) 
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and the perturbed stress-energy tensor for baryons contributes by 



ST° = - Pb 6 b , (39) 
ik l T° = ^ Pl 9 b . (40) 



All these ingredients are to be implemented in perturbed Einstein equations 



k 2 V - ~Kh = 4nGa 2 5T° , (41) 
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k 2 f] = AirGa 2 ik l 5T° , (42) 

h + 2Hh - 2k 2 T] = -%TiGa 2 8Tl , (43) 

1 

r 



h + 6rj + 2H(h + 61)) - 2k 2 r) = 24irGa 2 ( kjk,- - -5 i3 ) E) . (44) 



This system of differential equations can be integrated once fixed the appropriate initial conditions, which will be the 
content of the next sections. 

III. INITIAL CONDITIONS AND SUPERHORIZON EVOLUTION 

In order to start the numerical integration of the evolution equations wrote down in the previous section, one has to 
impose appropriate initial conditions to the fluid and metric perturbations. Although a general perturbation need not 
to be either isothermal (entropic) or adiabatic (isoentropic), it can always be expressed as a linear superposition of 
adiabatic and isothermal components p3[ . Also, it is useful to recall that isocurvature perturbations may be present 
in this kind of models [19] . We explore Doth these conditions in scalar field cosmology. In particular we search for the 
initial values of the fiela perturbations 6(f>o and 6<j)to (initial real time derivative) that realize initial adiabaticity and 
isocurvature, together with appropriate initial conditions on the other perturbed quantities. For this purpose, we will 
focus on the initial fluctuation of the real time derivative of the scalar field perturbation, since the conformal time 
derivative is always zero at a = by definition {5<f> — a ■ S(f> t ). The same happens of course to the initial conformal 
time derivative of the background scalar field 4>; if in general the latter has an initial non- vanishing kinetic energy, so 

that <f>to = (d(j)/dt)o 7^ 0, its conformal time velocity (j>o is zero since <f> = a ■ <f>t- 

In the following, we will need to use the scale factor behaviour at early times, when a« 1. We will often use the 
expansion of the scale factor in powers of the dimcnsionlcss parameter e: 



8irGp c tt r r— 

r = H y/il r T = Ct , (45) 

where p c and Hq are the present critical density and Hubble parameter respectively, and £l r = f2 7 + f2„ is the total 
radiation density contribution at the present; indeed, as it can be easily verified, being the scale factor a <^ 1 at early 
times, we can neglect the scalar field contribution in equation (Q), that admits the simple solution 

«W^+7?^ 2 + 0( e 3 ); (46) 



4 n r 



besides, the expansion rate behaves as 

U = C 



2 



1 1 n m 1 / f2 m 1 o 

7 + 4a--i6i7r) e + oie) 



(47) 



Before to go on, it is worth to remind some general results concerning the synchronous gau ge b ehaviour of metric and 
density perturbations on superhorizon scales (we refer to the work of Ma and Bertschinger |24j , although they did not 
include the scalar field component ). We impose initial condition at an early time, deep m the radiation era, when 
photons and baryons are tightly coupled and can be considered as a single coupled fluid; due to the large Thomson 
scattering opacity, the I > 2 moments of the photon temperature brightness function ( p4] ) (in particular, the shear <r 7 ) 

er in kr . one can ignore 
M), ©,@ become: 



and the polarization brightness function (25) are driven to zero; similarly, to the lowest or 



the I > 3 moments of the neutrino temperature brightness function . Thus, the equations 
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\k 2 S- ( = , (48) 

ifc 2 (^-4(7,) = , (49) 

^-{29 v + h + 6r}) = . (50) 

When wc impose initial conditions, at e <C 1, to get starting values for numerical integration, all the k— modes are 
still outside the horizon, i.e. k <C aH = 1/r (the last equality holds in a radiation dominated universe). Our aim is to 
extract the analytical time-dependence of superhorizon-sized perturbations at early times, once the initial conditions 
are realized: thus we find the early time form of equations (pl|)-(^4]), fl48|-|50|) , ( O ) and we find their solutions in 
successive powers of kr. To set-up the growth of perturbations, we must assume that at least a single perturbation 
is nonzero at initial time , in order to generate all the others. 



5v + ~9 u + ~h = Q , 6 V - 

2 

<i 7 - — (26> 7 + h + 677) = , ix y - 



IV. ADIABATIC INITIAL CONDITIONS 

The first necessary step to impose adiabatic conditions is setting to zero the initial entropy perturbation; ultimately, 
the origin of this result is that there is initially a single curvature perturbation (generated we suppose by inflation) 
and all later perturbations are inherited from it. The entropy exchange between any two fluid species a and b is ruled 
by the gauge invariant quantity 

S a b = 7- 7- , 51 

1 + w a l + w b 

that must be set to zero initially J19| . The second request comes from setting to zero the first time derivative of S ab ; 
actually, S ab obeys the following differential equation: 

Sab = -kV ab - iHT ab , (52) 

where T a b is defined as 

r a h = — — r a - — — r fc , 53) 

1 + w a l + w b 

T a being the gauge-invariant amplitude of the entropy perturbation of the fluid species a. The quantity V ab = v a — v b 
is the gauge-invariant difference between the gauge-dependent velocity perturbations of the species a and b. In order 
to have adiabatic initial conditions, both these terms on the right hand side of equation (^2|) are initially set to zero. 
Thus, for each pair of fluid components, we impose 

Sab = S ab - . (54) 

In particular, applying (|^) to the scalar field and another component (that we leave unspecified and label with x), 
will relate the initial values of 5<j> and 6<j> to the other energy component and metric perturbations. The first condition 



in (54) gives: 



< k8<f>t + V'6(t> = (ft T -r— ■ (55) 
l + w x 



Posing S<f, x — a ■ S^, x t = and using the Klein Gordon equation (|5|,pT|), we obtain 



1 



fc 2 t \ 1 



6a 2 HV J 6H 
Combining them together we find: 



1 M t - ( * + ( + 2V fc20? 



2 \l + w x J t l + w x \ a 2 V 



(56) 



(57) 
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6H - k 2 (j) t /(a 2 V') 





( 5x ) 







k + 



1 + 10, 



k 2 rh 2 

6H4h + 2V'-^ 



(58) 



The above expressions specify the general adiabatic conditions for the scalar field. Now, let us make a link to previous 
works; in ]24|] the adiabatic initial values and early time behaviours of the matter and the radiation components were 
found in the synchronous gauge; these results apply here, too. Indeed, as it can be easily seen from the Einstein 
equations, the contribution of the scalar field fluctuations is negligible at early times a <C 1 with respect to the matter 
and radiation ones, by a factor a 3 and a 4 respectively. Thus the approximations and treatment developed in ]24f| 
is valid also here for what concerns the ordinary fluid components, i.e. photons, massless neutrinos, baryons and 
dark matter; the time-dependence of the resulting super horizon-sized perturbations (fcr <C 1) is found by expanding 
the Einstein equations into powers of fcr and resolving the system of coupled differential equations to obtain the 
leading-order terms : 



8~j = 5 U = —6b = -6 C = — —Af{ki 



(59) 



15 + AR V 
23 + AR V 



<7 7 = , OV 



1 

~8 



45 + 12R V 



(kr) 



(60) 
(61) 



h=Af{kr) 2 , T) = 2Af- 



5 + 4R„ 
90 + 24i?, 



■AA(fcr) 2 



(62) 



where R v — p v /(p v + p 7 ) and Af is a normalization constant. Using these results, it is immediate to see from 
(§7|, !§ that, imposing adiabatic initial conditions, the initial values of 6<fi and 6<pt must be set to zero. Adiabatic 
conditions can be strictly verified only at this initial time, due to the effect of the mutual coupling between total 
density perturbations and entropy perturbations which appear in a generic multi-component fluid. Starting from 

initial zero values and using d62|), 5(f) and 5(j> will evolve according to the (O) which can be easily integrated once 
dropped terms of highest order m r; this gives the following behaviours at early times (a 1): 



^ t0 N Ck 2 r 4 



(63) 



1 



(f> t0 AfCk 2 T 3 



(64) 



having considered the lowest terms in r, thereby approximating the time derivative of 4> with its value at the initial 
time <fito ■ We have inserted these inputs into the standard CMB code and in Section VII we shall expose some 
numerical results. Now, let us turn to the second class of initial conditions. 

V. ISOCURVATURE INITIAL CONDITIONS 

The isocurvature initial conditions are obtained by setting to zero the gau ge invariant curvature perturbation. Its 
expression is given in terms of the gauge invariant perturbation potential \& plf| : 



2 (n- 1 ^ + ^ 

3\ 1 + w 



(65) 



we Point out that here, as in p9 24 1 , ^ indicates the gauge invariant <&a of the original work of Bardeen pH} |, while 
in |2lJ | the same quantity it is inalcated as $. Its expression in terms of the metric perturbations h and - ?} in the 
synchronous gauge is: 



1 

2fc2 



6rj ^ — ( h + 6fj 



(66) 
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Therefore, the appropriate isocurvature initial conditions are realized by the time growing solutions of the system 
©-© in which * and T^ 1 * are zero initially. First, let us see that, if the variables describing all the perturbations 
are regular enough to be derivable at least four times in r = 0, then the isocurvature initial conditions are simply 
imposed by setting the metric and radiation perturbations to zero initially. 

Isocurvature initial conditions : ho = r]o = <5 7 = # 7 — <r 7 = f>„ = U = a v = . (67) 

This can be easily seen by using essentially the Einstein equation multiplying both the members by a 4 , deriving 
once and factoring out the present critical density p c , it takes the Tollowing form: 

a 2 ^- (h + 677) + 4aa (h + 677) + (2a 2 + 2aa) (h + 677) - Ak 2 aari - 2fc 2 a 2 r) = -32irGp c {Sl 1 a 1 + Sl v & v ) . (68) 
Since by hypothesis h and r\ are derivable four times in r = 0, h + 677 admits the following early time expansion 

(h + 6r?)(r) = ^(h + 6r?) T + \^(h + 6r,) r 2 + ^(h + 6 V ) r 3 + ^—{h + 6r?) T 4 + 0(r 5 ) , (69) 



since with the initial condition ( |67j ) its initial value is zero. At r = the only term that survives in (|68|) is a (h + 6r)) 
since a§ = 87rGp c (^ 7 + n u )/3. Then, by using equations ( |i^ , ^9] , ^o| ) one obtains: 



(h + 6r)) ^^ Pc (Sly + n„) = (/i + 677) = . (70) 
In the same way, by deriving again (68) one gains 

% + 677) = . (71) 



Instead, it may be easily see that d 3 (h + 6rj)/dT 3 can be different from zero; in fact, deriving ( |68| ) three times would 
bring d?(h + 6t])/cIt 3 cx k 2 a 2 r\ that may be different from zero by hypothesis (|67|). This means that, for r — > 0, 

/t- + 677 = 0(r 3 ); since 7i = 1/r to the lowest order, this is evidently enough to make ^0 = (W~ 1 ^')o = Co = 0, showing 
that the initial conditions (J67J) imply isocurvature. 

It is evident that the initial condition ( p7|) can be realized in several ways, depending on which matter component 
is initially perturbed, or, in other words, on which S x is initially different from zero. In the present case a further 
degree of freedom arises from the presence of the scalar field, and we will analyze separately two main situations: in 
the first case only one matter component (CDM or baryons) is initially perturbed; in the second case the initially 
perturbed component is only the scalar field. 



A. Isocurvature conditions from matter perturbations 

Let us consider first the case in which an initial density perturbation, with amplitude 5 c q, resides only on the CDM 
component. By integrating (|35|), one finds 

5c = 6 cQ - h . (72) 

By hypothesis, this is the only initially perturbed quantity. All the others must to be set to zero at r = 0. Let us 
search the early time behaviour of the perturbations. Since all the modes are outside the horizon at early times, we 
first neglect all the terms proportional to k in the Einstein equations; then we expand all the quantities in powers 
of e defined in (^) and we calculate the leading orders. In making this, we are assuming that all the perturbation 
quantities admit a Taylor expansion in r = of course. By making use of the above criteria and of equations ( (46| , [47| ) , 
the Einstein equation ([Il|) becomes 

= (0^7 + ftv&v + n c S c e + Sl b S b e + 0(e 2 )) , (73) 



8 



and it is immediate to gain the early time behaviour of h: 

From the arguments exposed at the beginning of this section, up to the order e 2 we have also 

V = ~h. (75) 

Let us come now to the fluid perturbation quantities. As it is evident from the fluid equations, the 9 and a quantities 
are of higher order in kr with respect to the purely metric perturbations h and rj. Therefore, their early time behaviour 
can be written as follows: 

4 2 

&y = S„ — g^> = -^^h , (76) 



, a v = 0(e 3 ) 



(77) 



The behaviour of h + 677 is interesting even if of high order in r since it is directly related to the gauge invariant 
curvature by (pJ^Mfl) and it can be obtained by solving equation (144) : 



6r; = T r 



(78) 



where we have defined 



36ft r + 24f7„ 



(79) 



Note that ft. + 677 cx r 3 , according to the isocurvature nature of the present case, as we exposed in the beginning of 
this section. Also, ( |78[ ) can be used to find the behaviour of ov, by using again (RJ). 

It remains to findthe early time behaviour of the scalar field perturbation o^7 This can be done by expanding 
5<fi in powers of e and looking at the perturbed Klein Gordon equation once the terms proportional to k 2 have been 
neglected. The inhomogeneous term is — ^(frh; h is of the order zero from ([74]), and = a ■ 4>t is at least of the order 
e; thereby, to the lowest order in e, equation (111]) is satisfied by 



24^ 0<5c0 frC 



^ 1 c 2 



(80) 



This completes the early time behaviour of all the perturbation quantities in this case of isocurvature initial conditions. 
All these relations can be easily generalized to the case in which the initial perturbed matter component is the baryonic 
one. In the next subsection we study the other interesting case: the initial perturbed fluid component is the scalar 
field itself. 



B. Isocurvature conditions from scalar field perturbations 



Let us suppose that at the initial time a — > the only non-zero perturbed quantity is the scalar field, in a manner 
such that the total gauge-invariant energy density contrast is zero, all the others perturbation s be ing zero; this means 
that at least one of the two quantities ^Oi^o must be different from zero initially; from (|lq), the corresponding 
expressions for 6 (fro and 6<p t o are 



1 



V'{cfr ) 



<5 0O - 2-f ) + n0o)V> 



'HQ 



<frtoQ<t>o 



(81) 



In order to have isocurvature, for the other quantities we impose again the initial conditions (p7|). The relevant 
difference with respect to the previous situation lies in the slower rise of the metric and fluid perturbations starting 
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from their initial zero values: they will grow according to (f4l|)-(f44|) and ( |48| - |50] ), the whole perturbation-growth- 
machinery being initially driven only by the 0(e 4 ) contribution of the scalar field through the perturbed Einstein 

equations, while the perturbed Klein Gordon equation starts its dynamics from the conditions 5<f>o ■ 5<j> — and 
generates the inhomogeneous term driving the evolution of h. From equation (f4l|), together with (p^-pO|) , it is easy 
to find the early-time behaviour of the metric perturbation h: 

h = -( 8 jBI±\ C V + 0(r 5 } (82) 



4 \ p c VL. 

Using the method applied in the previous sections, one finds the leading-order behaviours 

4„ 4„ 2 



and, from (Eih it can be seen that 



where 



S~z = S v = -^S c = ^5 b = oc t 4 (83) 

„,0 7 ,0f, (XT 5 , (7„(XT 6 , (84) 



h + 6 V = ^r 6 , (85) 
o 



k2QA r^V-^-^nr). (86) 



170O r + %Vl v V Pc^r J V 10 10 

From the above formulas we see that the perturbations regarding the metric and the ordinary fluid components rise 
very slowly; indeed we found a substantial failure of this model in providing a significant amount of perturbations. 
For this reason we will not consider this case in the numerical integrations of section VII. 

It is interesting to find the behaviour of scalar field perturbation at early times, that moves it away from its initial 
value 6(f>o; this contains corrections in (kr 2 ) together with a term proportional to r 4 , as can be verified by integration 
of @: 

5<j) = 54> + H {2) t 2 + <50 (3) t 3 + 54>Wt* + <50 (5) t 5 + <50 (6) r 6 + 0(t 7 ) , (87) 
where the expansion coefficients are given by 

1 1 \l r 



(5) _ 



5 

"84 71 



-CS^ 5 



1 Q„ 

15 a 



i fn 



■CScjj^ 



i /n 



84 V 



c 2 - e 



80 V 

8<p^-C 2 V 



C 2 -k< 



J/o 



(88) 



We considered the expansion up to the sixth order in r because, as we will see in the next section, going to the 
Newtonian gauge changes the last coefficient. 

In the next section we extend the results of section IV and V to the conformal Newtonian gauge. 
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VI. RESULTS IN THE CONFORMAL NEWTONIAN GAUGE 



As it is well known, the synchronous gauge is a coordinate system corresponding to observers at rest with respect to 
the collisionless matter component. These "Lagrangian coordinates" are defined by the rest frame of a set of preferred 
observers. More physical intuition can be achieved in the conformal Newtonian gauge, where the metric tensor is 
diagonal. Inside the horizon, the perturbation equations reduce to the standard non-relativistic Newtonian equations. 
In this section we write the results of section IV and V in the Newtonian gauge. 

The connection between two gauges is realized in general by performing a coordinate transformation relating the 
two frames. The link between the perturbations in the two gauges is expressed in the same coordinate point instead 
of the same spacetime point; this is why in most cases it is interesting to know the difference of the fluctuations in 
the two gauges Q. 

First we writeaown the relations between the genuine metric perturbed quantities. In the Newtonian gauge the 
perturbation to goo exists and it is represented by the potential the trace of gij is instead perturbed by $. Their 
relation with h and rj are 



2fc 2 



h + 6ij ^ — ( h + 637 

a \ 



(89) 



* = 2PS(* + 6 *H- (90) 

They can be easily expressed for e <C 1,/ct < 1 by substituting directly the expressions for h and 77 contained in 
sections IV and V. 

Now we concentrate on the transformations regarding fluids and scalar field. They are the contained in the stress 
energy tensor, that transforms as 

dx^ dx v 

w dxP dx a w y 1 

Using this and taking care to compare the perturbations in the same coordinate point, the relations between the 
quantities in the two gauges (synchronous and Newtonian labeled as s and N respectively) for each fluid are: 

6 S = S N -T?-, (92) 
P 

9 s = e N - k 2 T , (93) 

Ps =Pn -pT , (94) 

a s = a N , (95) 

where 

- = ^ <«> 

is the lapse between the synchronous and the Newtonian time coordinate. Regarding the scalar field, we compute the 
Newtonian gauge expression of the amplitude fluctuation 8<p by using the transformation 



4>T . (97) 



In the following subsections we write the behaviour of the fluid quantities in the e -C l,kr <SC 1 regime and in the 
Newtonian gauge, dropping the ./V subscript. 



A. Adiabaticity 



The leading orders for matter and radiation perturbations are 

40A/~ _ 4 _ 4 

15 + 4fl„ ~ v ~ 3 b ~ 3 



*7 = - 1c , 75 =Sv=o8b=oSc, (98) 
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— <•»> 

k 2 r 2 . (100) 



, 2 2 



45 + 12i?„ 

The scalar field perturbation amplitude is 

^=^i5rk- (101) 

Note that in this case the scalar field perturbations grow in time faster (oc r 2 ) than in the synchronous gauge (oc t 4 ); 
as we point out below, this is not a feature of the isocurvaturc initial conditions. 

B. Isocurvature from matter 

Matter and radiation behave as 



5, = 8 v = \5 b = - 2 -8 c0 ^MT-^r. (102) 



5c = 5 c0 - \ 6 co^Ct - 3^ , (103) 

8, = 6 V = 6 b = -L 6c0 ^k 2 Cr 2 + ^f,e c ^^f. (104) 
a v = 0(e 3 ) . (105) 



The scalar field amplitude is given by 



■^c0^CV + ^«M, (106) 



and shows that the gauge change does not touch the order of the leading power in r, although it modifies its numerical 
coefficient. 



C. Isocurvature from scalar field 



Matter and radiation behave as 



^=^=-H^) cV - 2 !^ (107) 



^7 v v 



9 b cx r 5 , a v cx t 6 . (109) 



The scalar field amplitude is given by 



8<j> = 8cj> Q + 8<j>^T 2 + S^t 3 + 8^t a + S^t 5 + (V 6 ) + J|<M>) r 6 + 0(t 7 ) ; (110) 

therefore, in this isocurvature case, the behaviour of S(f> is affected by the gauge change only at high orders in r, 
leaving the leading terms unperturbed. 

In the next section we will numerically solve the linear cosmological perturbation equations with the initial conditions 
sets developed in the IV and V sections. 
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VII. NUMERICAL INTEGRATIONS AND DISCUSSION 



We performed the numerical integration applying our considerations to a scalar field model based on ultra-light 
pseudo-Nambu-Goldstone bosons; the potential associated to this field has the form [ p5[ : 

V{<t>) =A/ 4 [c OS (0//) + l] (111) 

Our working point corresponds to the parameter choice / = 1.885 x 10 18 GeV and M — 1CP 3 eV; assuming an 
initial kinetic energy equal to the potential one, the starting values of the scalar field and its initial time derivative 
are obtained by requiring that the present contribution is £2^ = 0.6, fixing Hq = 70 km/sec/Mpc and £2{, — 0.05. 
Furthermore, we have taken the primordial power spectrum to be exactly scale-invariant. 

Even though the main cosmological consequences of this kind of scalar fields have been analyzed by many authors 
(see [p|q]), here we use the formulas developed in the previous sections to accurately compare the pure adiabatic and 
pure isocurvature regimes. Also we give particular emphasis on the behaviour of entropy and curvature perturbations, 
again comparing their evolution starting from isocurvature (from CDM) and adiabatic initial conditions; each case is 
compared with a pure CDM model with the same background parameters. 

First, let us consider the power spectra of the microwave background anisotropies, both temperature and polariza- 
tion. They are expressed by the expansion coefficients of the two point correlation function into Legendre polynomials 
(see e.g. p3[) and admit the following expression in terms of the quantities defined in the previous sections: 

CT = ^ J f\A Tl (k,r )\ 2 ,£7f = 4tt J f\A Pl (k,r )\ 2 . (112) 

The adiabatic case is shown in figure 1 . The presence of the scalar field (solid line) produces an increase of the power 
of the acoustic oscillations with respect to the CDM model (dashed line) ; this is due to the fact that the universe is not 
completely matter dominated at decoupling: thus at this time the perturbations are growing faster than in the CDM 
models (we recall that density perturbations in adiabatic models grow as a 2 and a respectively in the radiation in the 
matter eras) and this produces an early integrated Sachs- Wolfe effect (ISW), found first in Q. Also, the position of 
the first peak is slightly shifted toward smaller angular scales due to the increase of the distance of the last scattering 
surface (projection effect). Note how these features regard both the polarization and temperature peaks. Finally, the 
temperature spectra shows that the ISW is active on the smallest multipoles due to the dynamics of the scalar field 
in the present case; this is a distinctive feature with respect to the cosmological constant 0. 

Figure 2 shows the spectrum from isocurvature perturbations. While the projection effect is the same as in the 
adiabatic case, now the situation regarding the amplitude of the acoustic oscillations is inverted: the peaks are 
lower than the ordinary models, both for polarization and temperature. This is simply due to the reduction of the 
matter/radiation ratio as we include <f> by keeping Qtotai = 1; in fact, the scalar field has no intrinsic dynamical 
effect at last scattering since matter and radiation components were largely dominant: it is well known the opposite 
behavior of the anisotropies in adiabatic and isocurvature models as one varies £2 m /i 2 (see e.g. Q). To better see this 
point, we plot in figure 3 and 4 the power spectra for models having fixed fl^h 2 and fl c h 2 but varying £2 m = fij + £2 C 
and h by mean of different £2^. Thus we expect the same amount of perturbations in the CMB except for the effects 
that are genuinely linked to the scalar field, as the projection effect and the ISW on the smallest multipoles. This is 
precisely what happens for the spectra in figure 3 and 4. The dashed lines represent again the curves for £2^ = .6 as 
in figures 1 and 2; the solid and thin lines represents £2^ = .5 and .7 respectively. Again, the spectra show remarkably 
the same features for polarization and temperature, even if it should be noted how the former, arising from acoustic 
oscillations occurring just at decoupling, is not influenced by the ISW effect, since it comes from the line of sight 
integration. All this shows how the perturbation power at decoupling is not touched by the subdominant scalar field; 
the opposite behaviours in the adiabatic and isocurvature cases is explained by the decrease of matter in favor of 
scalar field. 

More insight on the perturbation behaviours may be obtained by looking at the time evolution of some significant 
quantities; we look at one single scale, or wavenumber, roughly entering in the horizon between matter and radiation 
equality and decoupling: 

k = 8 x lO^Mpc" 1 . (113) 

Let us begin with the gravitational potential 'J, in figure 5; the oscillatory dynamics is associated, both in the 
isocurvature and adiabatic case, to the horizon crossing of the scale examined. The amplitude of the oscillations in 
the scalar field models are higher for the adiabatic case and lower for the isocurvature one when compared with the 
corresponding CDM models. These oscillations are the source of the CMB anisotropies on sub-horizon angular scales 
(I > 200) through the acoustic driving effect and the early ISW effect |l8j , and therefore follow the different behaviour 
in the two cases. 



We concentrate now on two particularly significant quantities regarding bo 
the total entropy perturbation, defined below, and the curvature £ defined in ( 
gauge invariant. The amplitude of the total entropy perturbation is given by 



adiabatic and isocurvature regimes, 
65J) ; we recall that these quantities are 
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pT = pT int + pT re i = ^2{Sp a ~ C 2 Jp a ) + X]( C a ~ C l) S P» = 
a a 

= Y J {5p a -c 2 s S Pa ) , (114) 



where p is the total pressure and the sound speed must take into account the scalar field contribution = 1 — 2V' 
in the summation: 



:> _ S a h a c a 



< = ' (115) 

where 

K = Pa + Pa , h = ^h a ■ (116) 

a 

While the Ti nt term comes from the intrinsic entropy perturbation of each component, and it is ultimately sourced 
only by the scalar field component due to the spatial and temporal variations of the scalar field equation of state, the 
T re i term arises from the different dynamical behaviours of the components, and it is related to the S a b quantities 
defined in eq. ([si]) by the relation Q : 

pr rei = i£^( c 2 -c*)s a6 (117) 

a, b 

In fig. 6 and 7 we plot r vs a (solid line and arbitrary units) using isocurvature CDM and adiabatic initial conditions 
respectively. We note that in both cases T — > as a — + 0: in the last case, this is an obvious consequence of the 
definition of what adiabatic conditions are; on the other hand, taking isocurvature CDM initial conditions, we started 
with non-zero values of the S a b relative to the CDM component and the other non-compensating components, but 
again r — ► as a — > 0. This is because at early times the cosmic fluid is radiation dominated, so that p oc 1/a 4 in 
( |l 14| ) ; this kills T for a — > 0, since no radiation perturbations are present initially. Instead, the initial value of the first 
time derivative of T is different from zero only selecting isocurvature initial conditions, since it takes contributions 
directly from the S a b terms |19[. The behaviours of entropy perturbations in both cases have been compared with 
those in the standard Einstein-De Sitter model fi m = 1 (dotted line), with the same choice of fib and Hq. The entropy 
perturbations remain nearly constant before horizon crossing; at this time the perturbation starts its oscillations that 
are damped in amplitude when the scale is well below the effective horizon. As an expected feature, note that in the 
scalar-field model the peaks of the oscillations are shifted closer to the present when compared to the Einstein-De 
Sitter case, as the epoch of matter-radiation equality. 

In fig. 8 and 9 we plot the evolution of the gauge- invariant curvature perturbation £ for isocurvature CDM and 
adiabatic initial conditions, respectively. At a <C 1, this quantity is zero in the isocurvature case and non-zero in the 
adiabatic one, being an explicit indicator of the nature of the perturbations. Again the significant dynamics occurs 
in correspondence of the horizon crossing, and the latter occurs slightly later then in the CDM model due to the 
presence of the scalar field. 

Finally, note how in all the isocurvature cases (figures 6 and 8) the amplitude of the oscillations are lower than 
in the corresponding CDM models; as we mentioned above, this is due to the reduction of the matter component in 
favor of the scalar field energy density. Most important, these graphs show that this is the only possible cause of 
this effect, since at the times of oscillations for the scale examined, roughly between equivalence and decoupling, the 
scalar field is very subdominant with respect to the other components. 

The hypothesis of a cosmic vacuum energy stored in the potential of a scalar field enlarges naturally the possibility 
to gain insight into high energy physics from the traces left in the cosmic radiation and in the matter distribution. 
Due to the upcoming CMB experiments [ p"7| , it will be interesting to further study the cosmological imprints of these 
models, in the context of different theories attempting to describe the hidden sector of high energy physics. 
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FIG. 1. Power spectra of the CMB anisotropies from adiabatic initial conditions. The background parameters are fib = -05, 
h = .7, three massless neutrino families and flj, = .6,fl c = .35 (solid line), Q$ = 0, fl c = .95 (dashed dotted line). Note the 
increase of the acoustic peaks power in the scalar field model. 
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FIG. 2. Power spectra of the CMB anisotropics from isocurvature initial conditions. The background parameters are Qb = -05, 
h — .7, three massless neutrino families and fi^ = .6, Q c = .35 (solid line), Q$ = 0, f2 c = .95 (dashed dotted line). Note the 
decrease of the acoustic peaks power in the scalar field model, an opposite behaviour with respect to the adiabatic case. 
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FIG. 3. Power spectra of the CMB anisotropies from adiabatic initial conditions with different Q$ and fixed Qbh 2 = .0245 
and n c h 2 = .1715 (fl m h 2 = 0.196) as in figure 1 and 2. The background parameters are Q. m = .4, h — .7 (dashed line), 
Q m = -3, h = .81 (thin line), and fl m = .5, h = .63 (thick line). The amplitude of the peaks is the same, while they are slightly 
shifted because of the projection effect. 
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FIG. 5. Gravitational potential (in arbitrary units) at the comoving wavenumber k = 8 x 1CT 2 Mpc -1 as a function of the 
time for adiabatic (top) and isocurvature (bottom) initial conditions. The background parameters are the same as in figures 1 
and 2. The oscillatory dynamics is associated with the horizon crossing of the scale considered. In the adiabatic (isocurvature) 
scalar field models, the oscillation amplitudes are larger (lower) than in the corresponding CDM cases, according to the power 
spectra behaviours. 



20 




FIG. 6. Gauge invariant entropy behaviour (in arbitrary units) as a function of the time in adiabatic models for scalar field 
(solid line) and pure CDM (dashed line) models. Note the shift of the horizon crossing (corresponding to the oscillations) 
toward late times due to the effective cosmological constant. 
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FIG. 7. Gauge invariant entropy behaviour (in arbitrary units) as a function of the time in isocurvature models for scalar 
field (solid line) and pure CDM (dashed line) models. Note the decrease of the oscillation amplitudes in scalar field models, 
due to the lack of matter with respect to the pure CDM case. 
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FIG. 8. Gauge invariant curvature behaviour (in arbitrary units) as a function of the time in adiabatic models for scalar field 
(solid line) and pure CDM (dashed line) models. Note that the curvature is non-vanishing aso^O. 
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FIG. 9. Gauge invariant curvature behaviour (in arbitrary units) as a function of the time in isocurvature models for scalar 
field (solid line) and pure CDM (dashed line) models. Note that the curvature is vanishing as a — » 0. 
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